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Converging:  Cylindrical   Shocks    In  Magnetohydrodynamics 

by 

G.    B.   V/hitham 


1.   Introduction 

In  this  report  the  method  developed  by  R.  F.  Chisnell 
(Ref  1)  for  treating  converging  cylindrical  shocks  in  ordinary 
gas  dynamics  is  applied  to  the  corresponding  problem  in  magneto- 
hydrodynamics.   Chisnell' s  method  is  approximate  and  relatively 
simple;  yet  it  leads  to  very  accurate  results.   For  strong 
shocks,  Chisnell' 3  predictions  may  be  compared  with  the  exact 
solution  first  calculated  by  Guderley  (Ref  2).   In  both  cases 
the  law  of  propa'-atlon  of  the  shoe'-'  is  found  to  be  U   r~ 
where  U  is  the  shock  velocity  at  a  distance  r  from  the  axis  of 
symmetry  and  n  is  an  exponent  depending  on  the  ratio  of  specific 
heats  y   ^^^  Chisnell 's  theory  n  is  given  by 


X. 


,i./ilpll,  (Z./^ 


The  values  of  n  in  the  tvjo  theories  are  shown  below  for  three 
values  of  y  (the  Guderley  values  are  taken  from  recent  caloula- 
tions  by  Butler  (Ref  3)  ) '■ 


Y 
n-  (Chisnell) 
n   (Guderley) 


1.2 

0.163108 

0.161220 


0.197070 
0.19729ij. 


1 . 666 . . . 
0.226220 
0.226051^ 


The  accuracy  of  the  approximate  method  is  well  confirmed 
by  these  results.   It  should  be  remarked,  however,  tb,9t  the  Gu- 
derley solution  gives  also  the  whole  flow  behind  the  shock 
whereas  Chisnell 's  theory  treats  only  the  immediate  neighborhood 
of  the  shock.   On  the  other  hand,  Chisnell 's  method  is  fiot  li- 
mited to  strong  shocks;  it  applies  over  the  whole  course  of  the 
shock  towards  the  centre.   This  is  an,  important  extension  since 
it  is  found  in  practice  that  the  final  period  in  which  the 


strongshock  assamptions  apply  is  often  quite  small.     The  pre- 
dictions of  the  velocity  and  strength  before  the  shock  becomes 
strong  have  been  checked  with  the  mamerical  computations  of 
R.  B.  Payne  (Hef  I4.)  and  good  agreement  kibs    found. 

In  the  corresponding  treatment  of  the  magnetohydrodynamic 
case  which  is  given  here,  a  further  advantage  of  the  Chisnell 
method  is  the  possibility  of  including  a  general  non-uniforrri 
state  ahead  of  the  shock.   In  the  Guderley  type  of  similarity 
solution  such  a  non-uniform  state  must  be  a  power  law  di.^tribu- 
tion  of  density,  pressure  and  magnetic  field  which  is  perhaps 
rather  artificial. 

The  Chisnell  rrtethod  in  general  deals  with  tre  propagation 
of  a  shock  down  a  channel  whose  breadth  varies  slowly  with  dis- 
tance along  the  channel.   The  converging  cylindrical  shock  pro- 
ble.-n  is  the  special  case  in  which  the  channel  breadth  is  propor- 
tional to  the  distance  along  the  channel.   i^or,  any  radial  line 
trirourh  the  axis  of  symmetry  is  a  pirticlo  path  and  could  be  re- 
placed by  a  solid  wall;  thu.s ,  if  we  consider  the  part  of  the 
shock  moving  down  the  narrow  sector  between  any  two  such  radial 
lines,  its  motion  is  identical  viith  that  of  a  shock  in  a  narrow 
wedge  shaped  channel  with  breadth  proportional  tc  r. 

The  propagation  of  a  shock  in  a  noxi-uniform  channel  is 
solved  in  the  followin-?  way.   Consider  first  a  channel  which  is 
of  constant  breadth  B  up  to  a  certain  point  but  then  changes 
its  breadth  by  a  small  amount  5B  after  which  it  is  again  ^oni- 
form.   If  a  shock  specified  by  the  pressure  ration  P  travels  up 
to  this  change  in  section,  a  transmitted  shock  with  modified 
stren.^th  will  continue  dox:;n  the  channel  and  a  weak  sound  wave 
will  be  reflected,   oy  an  appropriate  matching  of  the  various 
flow  quantities  (velocity,  pressure,  density,  etc.)  across 
these  waves  and  the  change  in  breadth,  the  change  in  shock 
strength  5P  can  be  expressed  in  terms  of  6B  in  the  form 

(1)  ^  =  f(P)6F  . 


Chlsnell  treats  a  shoi'jly  varyinii:  chaanel  as  a  series  of  succes- 
sive inf initessiraal  chanpes  5B  and  the  law  of  variation  of  P 
with  B  is  found,  therefore,  by  integrating  (1).   This  procedure 
neglects  interactions  bet'-jeen  the  reflected,  waves  from  succes- 
sive steps,  but  an  estimate  of  these  shows  that  their  effect  is 
small.   The  close  apreement  with  the  Juderley  solution  is  con- 
vincinp;  evidence  of  tiiis .   For  the  cylindrical  shock  B  is  pro- 
portional to  r  and  the  integral  of  (1)  gives  P  as  a  function  of 
r.   i'rom  the  shock  relations,  the  values  of  the  shock  velocity, 
density  ratio,  etc.  can  be  determined  in  terms  of  P. 

If  we  wish  to  include  a  variation  of  density,  pressure, 
and  man.netic  field  in  the  medium  ahead  of  the  shock,  jumps  in 
these  quantities  must  be  Introduced  wit?:  the  jiunp  /\^B  in  breadth. 
These  additional  chanpes  contribute  further  terms  in  (1),  but 
replacing  B  by  r  and  assuminf;  that  the  equilibrium  distribution 
of  density,  etc.  is  given  in  terms  of  r,  we  still  have  a  first 
order  differential  equation  for  P  as  a  function  of  r.   I''or  the 
limiting  cases  of  weak  and  stron^j;  shocks  this  equation  is  solved 
explicitly;  for  intermediate  cases  it  would  probably  be  necessa- 
ry to  compute  the  solution  numerically. 

The  investigation  given  here  will  be  based  on  the  Lund- 
qulst  equations  of  maeinetohydrodynamics  ('^ef  5)  ^nd  in  particu- 
lar tne  shock  relations  etc.  are  all  taken  from  the  recent  re- 
port by  Friedrichs  (Ref  6). 

2 .   Equations  of  iTiotion  and  the  shock  relations 

In  terms  of  the  velocity  vector  u,  the  magnetic  field  vec- 
tor H,  pressure  p,  density  p,  and  entropy  S,  the  Lundqulst  equa- 
tions of  marnetohydrodynamics  are 

(3)  li  +  "7  '^  (H  X  a)  =  0   , 

ik)  -^  +  (u-'7)u  +  i  T7p  +  txH  X  (V  X  H)  =  0   , 

(5)  It  •■  "7  •  (Pii)  =  0  , 


•.••■'■•     ;    I 


(6)  ||  +  (u  .  V)S  =  0   , 

(7)  P  =  P(P,3)   , 

where  ^i,  is  the  permeability.   The  Maxwell  equation  "V'H  =  C)  is 
essentially  included  in  (3),  since  the  divergence  of  that  equa- 
tion pives  '■H\7«H)/H  =0.   In  this  work  the  poly  tropic  law 

(8)  P  =  P^  ^(S) 

is  taken  for  the  equation  of  state  (7). 

For  flows  with  cylindrical  syinrnotry  we  assume  that  u  is 
a  radial  velocity  and  that  all  the  flow  quantities  are  functions 
only  of  time  t  and  distance  r  from  the  axis  of  symmetry.   -/'e  do 
not  exclude  transverse  and  axial  components  Hq  and  H   of  H. 
Indeed  thejr  are  of  primary  importance  since  it  is  easily  shown 
that  the  radial  component  H  ^  must  vanish  or  trie  solution  be- 
comes trivial.   F'or,  if  H   is  non-zero,  the  condition  that  the 
Q  and  z  components  of  H  x  (\7  k  H)  must  vanish  in  ( [| )  requires 
that  rKg  and  H   i.e  independent  of  r.   In  turn  these  lead  to  ar- 
tificial forms  for  u  and  so  the  only  feasible  case  turns  out  to 
be  Hq  =  H  =0.   But  then  the  flow  problem,  is  independent  of 
the  magnetic  field.   Thus,  we  consider  here  only  the  case  H  =  0. 

With  the  above  assumptions  of  symmetry,  equations  (3)  to 
(  S  .)  may  be  written 

(9)  -l+lfe  (Hq-)  -   0   ' 

(10)  -^-7-h  (-V)  -  ^    ' 

(11)    ^u  ^   ;>u  ^  1  ^  ^  1  it  ^  (^2  +  H^)  +  ii  -^  =  0   , 

^   '    ^-t  c;r   p  ,??r   2  p  -r  ^9    z'    p  r       ' 


^, 


(12)  |fi.,^.  p|a,£]i=  0 

(13)  If  H-  u||=  0  . 


;.ic 


If  p  (r),  p  (r),  H  (r),  refer  to  the  rvosslble  equilibrium  states 
that  may  exist  in  the  uadisturbed  medium  ahead  of  the  shock, 
these  functions  satisfy  the  relotions 

(1[|.)      p   +  p  li  H  +  ti,  1  — -—  dr  =  constant  , 

O     ^       O        v/    J^ 

(15)  S  =  constant    ,    Pq '"  Pq 

In  addition,  it  will  be  necessary  in  this  work  to  have  the  rela- 
tions connecting  u(r),  p(r),  p{r),  H(r)  in  steady  isentropic 
flow.   From  equations  (9)  to  (13)  we  may  write  these  in  tVje  dif- 
ferential form  (which  is  convenient  for  later  purposes): 

Hg  u 

H  u  r  ' 

z 

u  du   +   ^   +  ^    (H^,   dH^   +   H      dH    )    +   ii  H?,  —  =    0    , 
p        p'Q        Q  z        z  p©r  ' 

dp   +  du   ^   dr   ^    Q 
pur  ' 

dp  =  a  dp 
where  a  is  the  sound  speed.   Solving  these  equations  we  have 

/   6  _    du  _  -j^  ar 


Hg      u      r 

J  T  T 

dp  _  dp  _    z  _  ,,  , .  dr 

727  -  ~  -  "IT  -  ^^"^^  ■" 
a  p    ^  z 


) 


(I6)''j  P 

/  a^  +  tidip  -  Hg)  /  p 

where   A  =  —5 5 5 

a^  +  ^{eI   +  H^)  /  p 

Turning  now  to  the  shock  relations,  the  general  form  piven 
by  Friedrichs  simplifies  in  this  case  since  there  is  no  normal 


o 


component  of  the  nia.iiotic  field.   If  the  subscript  1  denotes  va- 
lues behind  the  shock  and  subscript  o  refers  to  the  equilibrium 
state  ahead  of  the  shock,  the  appropriate  relations  are 

p^(U  -  u-^_)  ==  PqU   , 
H-^dJ  -  u-l)  =  HqU   , 

12  12 

?!  +  ^  laH^  +  p^(u-j_-U)u^   ^  Po  "^  2  "^^^o 

1  ,,      ,2  ,     YPi     ^'4        1   2  ^     ^Po   ^  ^^o 

It  is  convenient  to  express  p,  ,  p^,  H,  ,  u-,  ,  U  in  terms  of  the 
state  ahead  of  the  shock  and  a  single  parameter  characterizing 
the  shock  stren^^th.   V/e  choose  E   =    p-j/p  as  the  parameter;  then, 

R^i-pf^TT?f%^i^o((2-y).  .r|] 

.  2Po^^--^^     r  2  ,  Y-1  >,2  ,^   ,s2, 
Pi  =  Po  -^  T^nr-T^-:!)!.   ^%  ^  V  ^o  (^-1)  ^  ' 

where  a   is  the  sound  speed  s/~Vp  /P   ^^^^   t)   is  the  Alfven 


speed  ■/[ill'^/p^      . 

In  the  limitinf?  cases  of  stron.^!:  and  weak  shocks,  the  re- 
lations (17)  take  simplified  forms.   Strong  shocks  (i.e.  large 

values  of  p-, /p  )  correspond  either  to  small  values  of 

10  2   2 

£  -  (y+1)/(y-1)  or  to  lar,"e  values  of  b  /a  .   The  former  is  the 

usual  case  found  in  jjas  dynamics  but  it  is  an  interesting  point 

that  in  magnetohydrodynamics  strong  shocks  also  arise  when  the 

magnetic  field  is  very  strong.  Ln   the  latter  case,  the  only 

simplification  is  that  a  may  be  set  equal  to  ?,ero  in  equations 

(17).   But,  when  £  -■>•  (y+1)/(y-1)j  the  shock  relations  take  the 

simple  form 


Y+1 


Pi    Y-1  ^O   '   -1    Y-1  -O 


H.  =  ^  H. 


(18)  { 

Pi  •   Y+1 


[ 


^o     r,2  2     ,, 

U     ,    U^  r:    -^  U 


It  should  lie  noted  that  in  this  limiting  case,  the  expressions 
for  P-,  ,  U-,  ,  p.,  in  terms  of  U  are  independent  of  the  ma(;riet.ic 


field. 


At  the  other  extreme  of  weak  shocks,  f  — >  1  and  we  have 
/ 


/ 


(19) 


/  p^  -  p^..-p^(?-l)   ,   H^  -  H^.^.  H^(^-l)   , 
u^,-  c^(r-i)   ,   p^  -  p^  ...■  Pq  a^(f-l)   , 


U  -  0   " 
o 


>^2 

(-H —  G   +  —r~ )  (£-1) 

'  ij.   o    Lj.   c^ 


where  c   =  /a  +  b 

^'      Interaction  of  a  shock  with  a  small  chan~e  in  channel 
breadth 

As  explained  in  Section  1,  the  main  step  in  the  Chisnell 
method  is  to  determine  the  strength,  etc.  of  the  shock  trans- 
mitted at  a  small  change  in  channel  breadth.   In  view  of  the 
subsequent  application  we  now  denote  the  channel  breadth  by  r 
and  consider  an  increase  5r  at  a  certain  section.   If  the  ini- 
tial  equilibrium  state  is  non-uniform,  the  ch-anfae  5r  i  s  accom- 
panied by  increments  in  p  ,  p  ,  H  .   Prom  (ll.|.)  these  increments 
are  related  by 


(20) 


/  6p   +  LiH  'SH   +  uH^^  -^  =  0 
'  o   ^— o  — o   "^  o9  r 


^o    o   ^o 


If  the  equilibriiom  field  H   is  civen,  p   and  p   are  determined 
^  — o  o      o 

and  the  increments  are  known  in  terms  of  5r. 

Let  the  incident  shock  be  specified  by  the  density  ratio  ^. 


!>■ 


r 

i.  - 


I'  '' 


8 

Let  the  values  of  flow  qvaantities  behind  it  be  u-j^ ,  p^,  p^,  H^, 
s3   =  YPi/Ptj  b^  =  ^H^/p-,;  these  quantities  are  all  determined 
in  terms  of  f  by  the  shock  relations  (17).   The  initial  situa- 
tion is  shown  in  Fig.  1.   ^fter  the  interaction,  suppose  the 
transmitted  shock  is  specified  by  the  density  ratio  £,   +    6^. 
Our  problem  is  to  determine  d^    in  terms  of  6r. 

The  values  of  all  the  flow  quantities  behind  the  trans- 
mitted shock  can  be  expressed  in  terms  of  £  +  5f:  and  the  equili- 
brium state  p   +  6p  ,  p  +  5p  ,  H^  +  5H^  from  the  shock  rela- 
tions (17),  and  we  denote  the  values  by  p^  +  Gp-j^,  p-j^  +  5p-|^, 
H   +  5H,  ,  etc.   There  is  also  a  contact  discontinuity  following 
the  transmitted  shock,  separating  particles  which  passed  through 
the  incident  and  transmitted  shocks,  respectively,  and  a  w:ak 
wave  is  reflected  back  down  the  channel.   In  all,  there  are  five 
different  regions  as  shown  in  Fig  2;  a  space-time  diagram  is 
also  included  in  the  firrure.   By  relating  the  values  of  flow 
quantities  in  adjacent  regions  in  an  appropriate  fasliion,  we 
shall  have  just  enough  conditions  to  determine  6F . 

All  the  quantities  in  region  o  are  known,  those  in  region 
1  (with  subscript  1)  are  known  in  terms  of  ^,  those  in  region  1| 
are  p,  +  5p,  ,  etc,  which  are  known  in  terms  of  r  and  5,^" .   V'e 
must  first  deduce  the  values  in  regions  2  and  3. 

Prom  1  to  2,  the  weak  wave  results  of  (19)  apply  ",  with 
a  m.odif ication  in  the  sign  of  u  since  it  is  now  measured 
in  the  direction  opposite  to  the  direction  of  propagation.   Thus, 
eliminating  £  in  (19),  we  have 

P2  -  Pi  =  -Pl(--~)  '   ^2  -  H,  =  -Hi(-^) 

(21)  '     2   u^^u, 

P2  -  Pi  "  -Pl^l  (-T^)   ' 

where  subscripts  refer  to  the  values  in  the  various  regions. 


Relations  (19)  actually  apply  to  any  small  disturbance  inde- 
pendent of  whether  it  is  a  compression  or  an  expansion  wave. 


9 


To  oass  from  region  2  to  3,  we  use  the  steady  isentropic 
flow  relations  (16),  so  that,  in  particular 

,     6r 

But,  u  13  continuous  across  the  contact  discontinuity  i.e. 
u-^  =  U-,  +  5U-,  .   Therefore, 

(22)  ^2  ""  ^1  "^  ^'^1   "^  '^l"!  'i^   * 

Now,  with  Up  and  u^  known  in  terms  of  f  and  5£  it  is  a  simple 
matter  to  find  all  the  other  quantities.   Region  2  is  determined 
from  (21)  and  region  3  is  then  found  using  the  relations  given 
in  (16)  which  may  be  approximated  to  first  order  in  6r  as 

I         '^91  "   "   1  r 

^^2        PjL        H =  (A-L  -  1)  -~   ♦ 

I 
The  final  condition  to  be  satisfied  is  that  p  +  p  l-"-^^-   ^^ 
cohtinuous  across  the  contact  discontinuity  ";  this  determines 
6E, .   We  have 

(21;)     p^  +  fip,  +  2   M'H^  +  I^Hi'SKj^  "^  P3  "^  i  '^^^3 

Prom  (21),  (22),  (23), 

^    f^      1^        2  5r  2  ,^^1  ^  ,   ^^1  5rv  _^ 

p  =  P2  +  (\^-l)p^a^  —  =  Pi  -  Pl^l  ^TT   '    h^T^    ^ 


'1     •"  "1 


/ ,   T  N    2  5r 
(A^-l)p-^a^  — - 


The  only  other  condition  across  a  contact  discontinuity  is 
that  u  is  continuous;  this  was  used  in  obtaining  {22). 


>,    I 


» 


10 

=  H  ,  -  H  ,  (— i  +  -1-i  •^}  +  (X,-1)H  ,  -^  . 
zl    zlc-,      C-,   r'     1   'zlr 

Therefore,  substituting  in  (2li),  we  have 

(25)  5p,  +  ^H^.5H,  -^  p,c^5u^  .  (^'^^  +  ^xH^^  ^^-i)  ^^  =  0  . 

All  the  quantities  in  this  equation  are  functions  of  £  given  by 
(17).   Hence  it  determines  5?  in  terms  of  5r. 

h'      Converg;in/^  cylindrical  shocks 

The  laws  of  propagation  of  converging  cylindrical  shocks 
are  obtained  by  integrating  {2^)    i.e.  ^(r)  satisfies  the  first 
order  differential  equation  obtained  by  substituting  (17)  in 
the  equation 

2 

(26)  r  |p  (Pi  *  |,h2)  .  p^c^r  ^  .  ^A^  .  ^H^^  ^M  ^  '  ' 

The  general  equation  for  ?  is  complicated  and  no  particular  pur- 
pose is  served  by  writing  it  down  here,  since  it  x-jould  have  to 
be  integrated  numerically.  Vie    can,  however,  get  the  solution 
explicitly  in  the  limiting  cases  of  weak  and  strong  shocks. 

Weak  shocks 

For  a  weak  shock,  the  approximate  relations  (19)  m.s^y  be 
used  in  (26)  and  we  obtain 

^  h   ^'o%    (^-1))  -^  Po^o^  h   (<^o^^-l))  ^  Po^o  (^-1)  -   °- 
The  solution  of  this  is 

(27)  ^Pq^o  (^-1)^  =  constant. 


11 


The  left  hand  side  of  this  equation  is  proportional  to  the  ener- 
gy flux  of  the  wave  at  distance  r.   Thus  the  law  of  propagation 
states  that  the  eneri;-y  flux  remains  constant  and  w  ^.  see  that 
(27)  is  the  anpropriate  generalization  of  geometrical  acoustics. 
From  (19),  the  variations  in  shock  velocity  and  pressure  jump 
are  given  by 

131 

Pi  "  Pq    "2  '2  "2 
^  p   c^   r 


(28) 


Po      °    ° 


'v  ,.  V.2    1   3   1 

j  ^  "  *^o    ,Y+1.  ^  2-Y   o.  "2  '2  "2 

V 

strong  shocks  with  g  "^  (  y+1)/(  y-l ) 

In  the  lim.itinp  case  of  a  strong  shock,  relations  (I6) 
may  be  used  in  (26),  and  we  have 

r  ^-    (p   u2)    +    p   U    /^  r  4"  +    p ^ =    0    . 

dr    '^o      '         ^o     /  Y-1        dr        ^o  /fTy^iy 


v/~2- 


The    solution   of   this    equation   is 
(29)  U  ci  p-''  r"P 


where 


(30)  a  =  ^^■— -      ,      P  = 


Y 


The  corresponding  law  for  the  pressure  p,  is 

/  -)i  ^  .^   1-2  a   -2(3 

(31)  Pi    Po    r      ^        , 

The  exponent  P  has  exactly  the  same  value  as  Chisnell 
found  in  the  absence  of  a  magnetic  field.   Thus  the  only  effect 
of  the  magnetic  field  for  strong  shocks  of  this  type  is  the 
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possibility  of  introducing  an  equilibrium  density;  distribution 
p  (r)  maintained  by  the  field  H  ,  and  satisfying  the  equili- 
brium relations 

1         2  ^^^^Go  Y 

p      +    o-  llH      + dr  =    constant  j  P^'^    P'  • 

^o    2  "^  o     r  o     o 

The  exponent  1  -  2a  in  (31)  is  always  positive  (provided  tnat 

Y  >  1).   Therefore,  if  p   increases  towards  the  center  p,  will 

increase  faster  than  in  the  case  H  =  G  and  vice  versa. 

P      2 
Strong  shocks  corresponding  to  large  values  of  ^2  /a 

If  ^  -  (y+1)/(y-1)  ^-3  i^ot  small  but  nevertheless  the 

shock  is  strong  because  the  magnetic  field  is  large,  the  shock 

relations  (17)  show  that 

(32)  p^  +  I  IXH^C^  p^b^   ,   Uc<  b^   ,   U^.v  b^   ,    C-^r/  b^   , 

where  the  factors  of  proportionality  depend  on  E   alone.   In  ge- 
neral, the  relations  are  still  so  complicated  that  an  explicit 
solution  of  (26)  seems  unlikely.   However,  in  the  special  case 
in  which  p   and  b   are  powers  of  r,  it  may  be  observed  that  the 

•^O  O  7  O 

terms  in  (26)  are  homogeneous  in  r.  Hence  there  is  a  solution 
of  the  equation  in  which  f  takes  a  certain  constant  value.  In 
this  case  the  dependence  of  p,  ,  U,  etc.  on  r  is  given  by  {32.}. 
The  existence  of  this  type  of  solution  was  pointed  by  Cathleen 
S.  Morawetz  who  first  found  it  in  her  investigations  of  the 
similarity  solutions  of  the  full  equations. 
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